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I. INTRODUCTION 


This Is an InterrIm report of the second phase of research under 
Grant NSG 1128, dealing with the design of a Microwave Landing System 
(MLS) aircraft receiver capable of optimal performance In the multipath 
environments found In air terminal areas. The project focuses on the 
angle-tracking problem of the MLS receiver; the work reported here 
includes tracking system design considerations, continued study and 
application of locally optimum estimation. Involving multipath adaptive 
reception and then envelope processing, and finally microcomputer system 
design considerations. A significant result obtained is that envelope 
processing is competitive In this application with i-f signal processing 
performance-wise and is much simpler and cheaper, hence future effort 
will focus on envelope processing. To provide a basis for discussing the 
results obtained, a brief summary of the signal model concludes Ihls I 

Introductory section. 

In the last report [1] the received signal In the (linear) i-f 
channel was modeled as a function y(tj^,T) of a global discrete-time 
variable t^^ and a continuous time variable t local to the present scan, 
as fol lows: 

corresponding to direct-path and reflected components and receiver noise, 
respectively where, neglecting doppler effects, on the kth received scan 
and for 0 i t s T, 

y (t ,t) = a(t )p[6 (t) - 0(t. )]coi[(D x + 6(t, )] (I-?) 

uk k." k Xt* ^ 

YR^+k^T) = I Oj(+^,t)pC0^(t ) - 0p (t^^)]cos[o)jpT + (1-3) 

= y^ (t|^,T)cosCa)jpT + * Yp (t^,i )si nCu^pX + S(tj^)] (1-4) 

c ^s 



n(T) = stationary, bandpass, Gaussian process, mean zero, 

variance o ^ 
n 




= n (t)cosCu T + B(t )3 - n (T)sinC(*>_cT + 6(t, )3 (1-6) 

c IF k s IF k 


Alternatively, in terms of the above, we may write also 


y(t|^,T) = V(t|^,T)cos[u)jpT + 6(t|^) - r(t|^,T)] 


(1-7) 


where 


V(t.,T) = /{a(t. )pC0. (T)-0(t, )D+y„ (t ,t)+n (T)}^+{y„ (t ,r)+n (t)}^ vI-H) 
K KM k K K C k 5 


Yr (t(^,T)+n2(T) 

r(tk»t) arc (t|^,T )+n^(T ) ^ 

c 


(1-9) 


Ar,(t ) 

Yr ^ ^ “i ' c^^ ‘ ®R ^ + (1-10) 

c i i 


Ar. (t ) 

Vr [• “i ^■'■k'‘^^Pl-®A^‘^ '~c^^ ‘ ®R ^^k^ ” (I-ll) 

s i ^i 


6i(t,) - 6(t^) = e.(V,) - 6(t,.,) - f (t, - 


Refe ence to Cl3 is made for definitions of unfamiliar quantities not 
defined (or redefined) above. The principal parameter to be estimated is 
the aircraft (A/C) angular coordinate 0(t. ); other parameters whose estinates 

K 

are needed for the 0-estimation include the amplitude parameter ai'id 

the noise variance Necessary modeling of the evolutionary dynamics of 

0 and a has been deferred until the requirements of ftu^ osfimation a I 'K'r i I hm 
chosen become firm. 
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II. TRACKING ALGORITHM DEVELOPMENT 


The sought tracking algorithm for the MLS receiver is a discrete- 
time estimator of a vector x(k), comprising the A/C angular coordinate 
and ancillary variables, given the sequence of observations (V(k), k=l,2,...}, 
where, for the kth scan, 

V(k) ^ {y(tj^,x), 0 < T < T}. (II-l) 

Such algorithms generical y are characterized by the following functions: 
Extrapolation : xfk|k-l) = f (x(k-l |k- I ) , k, k-l) (II-2) 

Error Estimation : e(k|k) = g(V(k), VCklk-l)) (II- 3) 

Updating : x(k|k) = x(k|k-l) + e(k|k) (II-4) 


where 


x(k-l|k-l) - an estimate of x(k-l), given all observations up 

through the (k-l)th scan. (II-5) 

x(k|k-l) - the extrapolation of x(k-l!k-l) up to the beginning 

of the kth scan (II-6) 

e(k|k) = an estimate of the error in xCkjk-l), given the kth scan 
observations V(k) (II-7) 

V(k|k-I) = predict'on of V(k), based on the extrapolated estimate 
x(klk-l). (II-8) 

The dependency of e(k|k) on x(k|k-l) through V(k|k-I) indicates that such 

algorithms are recursive. The theory and design of recursive state esti- 

mators is we I I -documented, and given suitable, valid models , this approach 
might be applied to the extended problem involving also identification of 
several model parameters Imprecisely known. Modeling the state excitation 
as white noise and augmentation of the state with the parameters to be 
identified are required generally. Recursive estimation has much to 
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recommend It generally, but the potentially high dimensionality of the 
augmented state model and reservations about the validity In this appli- 
cation of modeling state evolution uncertainty with a white state noise 
both have made more attractive a layered approach to the extended problem, 
as fol lows: 

1. Use modified recursive estimation of the angular coordinate 
and ancillary variables (simple state). 

2. Use batch (l.e. finite-memory) processing of a sequence of most 
recent, raw state estimates to obtain a more refined state 
evol utlon model . 

3. Extrapolate with the most recent refinement of the state 
evolution model . 

Modified recursive estimation refers to the error estimate e being defined 

and calculated for its direct addition with constant unity gain in the 

update equation. This is necessary to insure that the results of the 
batch processing include the most recent observations, fully weighted. 

In this way also both tracking stability (including false lock) and esti- 
mation quality (including suppression of multipath effects) can be dealt 
with in developing the estimate e. For example, clearly e(k|k) is a 
function of the actual error e(k) in x(k|k-l), where 

e(k) ^ x(k) - x(k|k-i). (II-9) 

Strict stability of the tracking algorithm requires 

<e> = 0 when e = 0 (II-IO) 

0 < <e> < 2e when e > 0 (II- 1 1 ) 

2e < <e> < 0 when e < 0. (II-I?/ 

Estimation criteria for e should optimally relate to the quality of the 
X estimates, perhaps, but with the assumption that x is a fixed parameter 
throughout the kth scan observations ^(k), a reasonable, albeit possibly 
suboptimal, approach is to require e to be locaily unbiased at e = 0 and 
have minimum mean square error, i .e. 
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= 0 (vector) 


(II-I3) 


Local ly Unbiased : <e - e> 

e=0 



(matrix) 


(II- 14) 


MMSE: 



<(e - e) (e - 


e)^> 


e=0 


be minimum among all estimates (u_| 5 » 
of e locally unbiased at e=0 


This is the locally optimum estimation criterion expounded by Murphy in 
[2j aiid applied in ClD and in the next chapter of this report to estima- 
tion of the A/C angular coordinate, given respectively the (linear) i-f 
signal anc Then its envelope. Approximate covariance-of-error expressions 
for the estimate e are produced also in terms of the estimated signal-to- 
noise ratio on the kth scan; these might be useful for weighting the raw 

A 

es*‘imates x in the batch processing, as shown below. 

Batch processing, corresponding to finite-memory or moving-window 
filter !''q, is used here to fit in a least squares sense on assumed state 
evolution iv.'-* linear in the unknown parameters, constituting a vector V, 
to a sequence of estimates x, constituting a vector X, via minimization 
of the quadratic form 

(HV - X)%(HV - X)| = (WHV - WX)^(WHV - WX) (II-I6) 


giving 


V = (WH)\x = V(k) (II-I7) 

where ( )^ denotes the pseudoinverse of the matrix ( ). If a linear law 
of evolution is assumed of the form 


x(t) = x(t|^) + (t - t^)x(t^) 


(II-I8) 
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and the most recent K + I raw estimates x and associated error variances 
0 ^^ are used, then 


X = 


x(k |k) 
x(k-l |k-i) 

x(k-K|k-K) 


H = 




I <VK-+k> 


and 


k 


Tin 


w H 


a (k-l) 
e 




\ ' i / 

\ 0 • • • ^nr-"KT / 


(11-19) 


( 11 - 20 ) 


(II-2I ) 


Under these circumstances H and WH are K + i x 2 with full rank 2 and 
thus as an alternative to the expression (II-I7) for V above, we have 


V 


(h\h 




V(k) 


( 11 - 22 ) 


In all cases the matrices to be inverted are either 2 x / in dimensicr or 
otherwise diagonal (approximately), so the computational load for model 
estimation is not excessive. The number of measuromonts k(> ?) used here 
needs to be established; this might be done in a manner that would make 
the algorithm somewhat adaptable to A/C maneuvering. This will be studied. 


i 
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Given a solution for V(x), extrapolation for x(k+l |k) is nccomp I i 
as fol lows: 

x(k+l Ik) = (I (t. ,-t, ))V(k) (11-73) 

' k+l k 

Wo note this Isa smoothed prediction; the smtxjthed "current estim.jte" 
value is available as the first element of V(k). The one-step prediction 
Is exoected to be the principal result, however, both for maintenance of 
th • algorithm (In calculating e) and in output, since processing time is 
off.tt somewhat In prediction. 

In summary, a layered tracking alrorithm structure has been described 
invciving recursive estimation of the A/C coordinate and batch processing 
of these estimates for model identification. The approach allows false— 
lock prevention and suppression of local multipath effects (assymetrica I 
pulse distortion) to be Included in the well-defined problem of processing 
new observations for optimal estimation of prediction error. Also the 
approach permits the batch processing window to be adaptable to manifest 
A/C dynamics, thus producing a good and reasonably recent model for 
extrapolating the estimate through gross multipath effects, such as sianal 
fades. A disadvantage of batch processing is the storage reguirement for 
past data, if significant. An overall study of processing time, storage 
regulrements and tracking performance is being done. Ful I y-recursi ve 
approaches have not been totally discounted either and one also being 
studied. Particularly attractive are the two-filter versions of Biemian 
[3] for fixed memory filtering and Nelson and S1ear [4] for simultaneous 
state and parameter estimation. Batch processing and recursive estimation 
approaches ma, not produce the same results [53, and performance analyses 
and comparisons are essential, along with assessments of computational 
loads. In selecting the algorithm to be implemented for field test. 
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III. LXALLY OPTIMUM PREDICTION ERROR ESTIMATION 
Part A: BASIC THEORY AND l-F SIGNAL PROCESSING 


The concept and development of locally optimum esilmation, expounded 
by Murphy In 1968, was summarized in CO and applied to the A/C angular 
coordinate estimation problem In MLS, given the receiver i-t signal. To 
facilitate presentation and evaluation of new results to be given, .1 brief 
summary of some results from CO Is presented in tbi'. ch.ipter including 
key relations In the locally optimum estimation model. 


Locally optimum estimation theory involves operations on the Radon- 
Nikodym derivative dP /dQ where P is the measure corresponding to tfie 
integrated observations process (Y(t), t c C0,Tl), and 0 is the measure 
corresponding to the integrated receiver noise process {N(t), t e C0,T]}; P^ 
is absol utely-conti nuous with respect to (wrt) Q. The Radon-Ni kodym deriva- 
tive is a generalization of the ratio of the two approfjriate probability den- 
sities, commonly termed the likelihood ratio to which it degenerates when * 

tfiR measures P and Q are both absol ute I y-cont i nuous wrt I obesque moasurc p . 

0 

Without regard for these finer distinctions, the symbol X and term likeli- 
hood ratio will be used indlscrir-.inantly for either mathematical object, 
as applicable, thereby rendering the basic model applicable to both con- 
tinuous-time and discrete-time estimation, the latter being used in the 
envelope processing algorithm where probability densities of the (finite) 
sets of samples are available. Specifically, therefore 


, for contI nuous-t Ime processes with 


LIkol i hood Rati 


o. = I 


measures Pg»0 as defined. 

I ^ * prese n t ^ discrete-time processes, 


sigi.al absent 


V being the K- vector of 
envelope samples, , 

CV(t|^»T|), ... V(tj^,T|^)J . 


(IIl-l ) 


Murphy [2.1 has shown that the estimate e of the vector e which is 
locally optimum at e = 0 Is given by 


0 



. -1 
e = ♦qAo 


(I1I-2) 


whore Aq 


Is the vector whose 1th component, 



I s given by 


A 


0 


I 



0=0 


if A / 0 
o 

otherwise 


(II1-3) 


and 


^A()Aq> , 


(1II-4) 


the latter expectation being taken wrt the measure corresponding to 

e = 0. Further, If e = 0, the residual mean square error (I.e., the 

-1 

error covariance matrix associated with this estimate Is ♦q , that Is 


• aT I 

<ee > 

'e=0 


-1 

♦o . 


(Ill-*;) 


Applying these results to estimation of the error In thu one-step 
prediction x(k|k-l), given the kth scan observation l^(k), the error vector 
e was constituted as follows: 


e 




e 





error In the prediction, 0(t 
error In the prediction, a(t 


k 

k 




(HI-6) 


corresponding to the state vector x being estimated. Four cases have been 
or are under consideration at this time, as follows: 

1. I-F signal observations 

a. No mul 1 1 path 

b. With multipath 

2. Envelope observations 

a. No multipath 

b. With multipath 
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These studies and results are descril)od below and in the following chapter. 
I . I-F Signal Observations 


The i-f signal y(t. ,t) was given in Chapter I as follows: 

K 

y(tk,T) = + n(T) (III-7) 

where, for 0 < x < T on the kth scan, neglecting dopp I or effects, 

y_(t. ,t) = n(t. )p[fl.(t) - 0(t. )]cos[a»,_T + h(t, )] (IIl-R) 

D K K A k I r K 

' I “I'V-apCe^iT - * ®i‘V’ 


and, since the i-f bandwidth B(= 160 kHz) is suhstan t i a I I y greater tnan 
the reciprocal puise width (8 milllsec"^. 


n(t) = white Gaussian noise with 2-sided spectral density N (' o^/2B). 

on 

(III-IO) 

in the following y^ denotes the above expression for with the estimates 

a(t It ), e(t |t ,) substituted for a, 0 respectively; also the explicit 
k ' k- 1 k ' k- 1 

dependency on t^ is suppressed when no confusion results. 

a. No Multip ath. Hero yp = 0; under the hypothetical assumption 
that the pari».meter 6 Is known, the likelihood ratio. A, for instrumentation 
purposes, can be written as 


A 

e 


= exp{^ yp(x)[y(x) - >iy|^(T)]dx} 


(III-I I ) 


Also, under the same circumstances. 


Ac = ^ 

o 


yp(x)]dx 


/ T ^yn ^ ^ ^ 
0 60 
,T 


(III-I2) 
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where 0, a denote 0(tj^ respectively. Finally, and most 

Importantly, after substituting for y (III-8), and some manipulatior in.i 

-1 ^ 

simplification, we obtain ♦q» covariance of the estimation error 


c 


a - 

= e - e 


for e = 0, as fol lows: 


-1 


{ot^ p2[0 (T)HdT) ^ 

■' 0 A 


= ?.H 


O A.n 

o 


P^C0y^(T)]dT> ^ 


where 


pCe^(T)] = ^ pfQeH 

le = 0 „(t), 

e A 


(III-I3) 


(ITT-14) 


(III-I3) 


The case for B unknown but assumed to be a random variable uniformly dis- 
tributed on (-• iT,Tr) will not be summarized here; reference is made to [ I ] 
for this develooment and the resulting structure, which is a quite compli- 
cated quadrature detector signal processor. 

b. With Multipath . The processor of the signal for the no-multipath 
case above, as evidenced in the expressions for Ag above and in [l], consists 
of an integration over the scan interval of the product of the received 
i-f signal and a sinusoid which is both phase-synced with the i-f signal and 
amp I i tude-modu lated by the derivative (wrt the parameter of interest) of 
direct path envelope function. This Is a form of gating (weighting^, c'jr 
and in reaching for a concept of multipath-adaptive reception, it seemed 
reasonable to expect the same general structure in the nucleus of an adaptive 
receiver with the amplitude-modulation of the local oscillator being con- 
trolled by suitable additional algorithmic machinery to adaptively discrimi- 
nate against manifest multipath interference. These notions are fairly 
easily confirmed for the i-f signal observations case, using the locally 


optimum estimation equations. 

The multipath propagation component the composite I-f 

signal y(t|^,x) is modeled as the sum of Individual reflection components 
y (t ,t), each a function of an amplitude variable a.(t ,x) and a phase 

r\ j I h 

variable B.(t ). Within a given scan period the B. have been effectively 
IK I 

determined to be independent random variables uniformly distributed on 

(through their Individual variations from scan-to-scan are highly 
correlated, thus Insuring the continuity in time of the interference 
phenomenon present). Under the assumption each a.(t ,t) is a (time-varying) 

I K 

Rayleigh random variable (i.e., a Rayleigh random process) independent of 
the associated B.(t ), each reflection component y (t ,x) is a nonstationary 

IK ^ j K 

(jaussian process, as is consequently also the total multipath propagation 
component Under milder assumptions on the a.(t ,x), if the 

number of non-zero ’s are large, it is still possible to argue that 
y^(tj^,x) is approximately (Baussian. It seems plausible that usually some 
combination of the preceeding conditions would prevail, such that the 
multipath interference component yp(tj^,x) in the received i-f signal 
y(t^,t) is an additive (Gaussian random process with zero mean value, or 
more appropriately a random field with zero mean value and covariance 


This assumption we make. 

Let the composite i-f signal y(t ,t) be sampU^d on the kth scan 

K 

and a K-vector Y(tj^) be defined, comprising these samples, i.e.. 




1 





1 

Y(t, ) = 
k 

1 • 


• I 

+ 

VR(+k-''2> 

1 • 

+ 

n(x2) 1 

:• I 










= Yo(tk) - - N], 


(III-I7) 


(III-I8) 
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defining the K-vectors Y|^, Y^, and N. Regarding Y^ as the signal (known 

sure function) and Y + N as the interference (Gaussian), then the likoli- 

K 

hood ratio, corresponding to (III-I6) above, is 




p (Y) ' exp(-4[(Y-Y^)^(t^*»^r'(Y-Y^) - y'"(V«^)"y:|I (I1I-I9) 




where is the covariance of the noise vector N, and 4>„(t. ) is the co- 
N R k 


variance of the reflection component vector Y (t. ) 

R K 


(The ijth entry in 


♦„(t ) is ♦oft ,t ,T,,T.), from (III-I6).) The above can be simplified to 
R k K k k I j 


-1 


X = exp[Yp(4>^ + 4>|^) (Y - bY|^)], 


(III-20) 


corresponding to equation (III-II) above. Then, corresponding to (III- 1 2), 
the optimal processing of the received signal Y(tj^) is indicated by 


/ -T \ 

(*R + $.) (Y - Y ) 
30 R N D 


^0 


\ 


'■T 

-1 . 

— (4>„ + ♦„) (Y - Yr,) 

33 R N D 


(III-2I ) 


which, if the i-f noise samples are independent with variance o^, then 

n 

4>,, = 0 ^ 1 , and An can be expressed in the form 
N n “ 



[I 




\ 


da 


[I 



n 


(I + — 


-1 

) ] (Y 





R 


-1 


+ -T ) ] 


(Y 


\ 



(III-22) 
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The innovations (Y - Y^) are effectively processed wifn a welcihtlnii t imc t iv’'n 

3Yr\ ^ 4^ 

D , . R . R 


that is dependent on 1 .e. , -jr- I 


-1 


1 - -y (I + -j 


) ], In the aninil.ir 


coordinate channel. When there Is no multipath and hence 4*^^ = 0, the 

results correspond with equation (III- 1 2) above, oth^'rwlse a modification 

dependent on $ is made that preserves locally optimum estimation performance. 
K 

Application of the algorithm requires estimation of t , of course, as well 

R 

as t,,; however as (III-2I) shows, only the sum + 4>^,) is needed not 
N R N 

the individual covariances. The sum (4> + 4> ) Is the covariance of the 

R N 

Innovations process (Y - Y^) at e = 0, however, suggesting a method by 

which + 4>., might be determined and even tracked as 4»_ varies with the 
R N R 

environment. 


No further study of i-f signal processors in general, and their 
multipath suppression capabilities in particular, h.is been done, or I;, 
planned In view of the competitive performance observed of the much more 
economical envelope processors, but the notion of interference-adaptive 
reception is clearly capable of generalization to non-additive, non-Gaussian 
interference models. This work is In progress. 
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IV. LXALLY OPTIMUM PREDICTION ERROR ESTIMATION 
Part B: ENVELOPE SIGNAL PROCESSING AND PERFORMANCE EVALUATIONS 


As pointed out In the December 1975 interim report [l], the com- 
plexity of the optimal l-f signal processors motivates serious considera- 
tion of the simple envelope detector and the optimal processing of the 
resulting envelope samples In the estimation of A/C angular coordinates. 

To date the envelope processor has been developed only for the no multipath 
case and an outline of this development Is presented below. 


Murphy C23 (see also Chapter I 1 1 of Cl]) has shown that the locally 
optimum estimate e of the parameter vector e at e = 0 Is given by 


-1 

e = ♦qAq 


(IV-I ) 


where Aq Is the vector whose ith component Is 


Ao. =' 


dP dP^ 

I 'e=0 


, otherwise 


(IV-2) 


and 


- <AqAq> = 


/ AgAodPo 

n 


(IV-5) 


It is assumed that, for the samples of the envelope, the Radon- 
Nikodym derivative is equal to the likelihood rafio, i.e.. 


^*^0 _ p(v|s+N) 

e=0 " 


e=0 


(IV-4a) 


= Io(-^y-)exp(-^) (IV-4b) 

J=l *^n 2a2 

n 
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whicn further assumes the 1-f noise is a narrow h.ind, <'ero mean r>Rf' with 
variance and that the envelope detector function is c)ivon by t\iuafien 
(1-3) with y = y_ = C These latter assumptions lead to the detormi- 
nation that tne prooability density function for the envelope samples 
given that e = 0 is C p. 166 and p. 357] 


V. 


V.P. 


Vt + Pt 


p(V,le=0) = 4 io(-i-J.)exp - {.I - 4 ) 


2ai 


V > 0 
k ~ 


(IV-';) 


where {Vj} is a set of identically distributed (Rician) random variables . 

Note that V. = V(t ,t.) and P, - Pit^,T.) = a(tj^)p[0^(T . ) - 0(t^)3 = ap . 

J J J J J J 

are the sample values of the envelope and the deterministic signal due to 

the antenna selectivity function, respectively, at t = t. on the kth scan. 

Also recall Cll that e is the error vector in the parameters being estimated, 

and in the present case 


e = 


^0 


(IV-6) 


where e =0-0 and e = a - o, i.e., ■'"he errors in estimating the 
0 and the amplitude of the direct path signal a; and therefore e = 0 implies 
0=0 and a = a (not 0 = a = 0). It should also be noted at this point 
that must be assumed known (or being provided via a separate estimation 
process) . 

The form of the likelihood ratio given in (IV-4b) is obtained by 

assuming independence of the K envelope samples which implies p(Vj, ... V^^) 

K 

= IT p(v.) and dividing by the product of Rayleigh densities p(V.|P.=0K 
j=l J J J 

The independence assumption implies that the sampling interval At = i 

Is equal to the value of the delay variable which makes the autocorre • a''" i cn 
function of the i-f noise equal zero [7, p. 399 and p. 416]. If the i-f 

- 1 

noise spectra is assumed to be ideal bandpass, then At = (i-f Bandwidth) 

■ 6ys. 
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The components of Aq ^rom (IV-2) aro the partial derivatives of the* Kui 
likelihood ratio and are given by 


A k 
" 30 


e=0 



P. } 

J 0=0 


(IV-7) 



e=0 






P.} 

e=0 


(IV-8) 


(Note that qj = P 2 = 0 are the maximum likelihood equations for , 

I.e., 0 = 0... and « - u,., .) Thus 
ML ML 


Aq^o - 


Q Q 


q q 

2 1 




and 


♦o - <AqAo> = 




<q q > 
1 2 


<q2> 

2 


(IV-IO) 


As the locally optimum estimator e is unbiased at e = 0 


-1 -1 
<o> = 4> <Aq> = 


A 

n 

V 


[0 ' 

1 

- 

0 



2 j 

1 

\ 1 


nv-i I ) 


To show that this is indeed the case, it suffices to show, using frec-lv 
some notctional abuse. 
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v.p, 

Il(-y) 

i„(^) 


(iv-i:) 


= r V {_ 

^ j io< ) O 


V + P? 

{-4 Io( )exp - (J ^ 


‘■°n 


V2 + p2 


= r -4 Il( )exp - ( J ^ civ. 

« °n •' 


= P^.Q2(x,0) = P^. 


where 


A Z ' 

Q (x,y) = / 2(-) 
m ' ' y X 


I ,(xZ)exp 
fn— I 


Z2 + x2 
( ■ 2 ) dZ 


(IV-13) 


is the generalized Marcum Q function [8, p. 41 lU and 0 (x,0) = I, all m 

m 

and X (the Integrand is the probability density function of the random 
variable Zo = |1 y|| where Y c (A,o^) and xo = I|a 1|, see C*5, pp. 4l-4?3) 


where Y c N 2 ^(A,o^) and xo = 


Now the expressions for and (equations (lV-7) and (lV-8)) may 
be simplified by letting 


y ^ V llLl. p 

j Iq( J j 


(IV-14) 


be a new random variable with <y.> = 0 and <y,y.> = 0 for i ^ j. Thus 

J i J 

(lV-7) and (lV-8) may now bo written as 


' r • 

~ 2. “PiY; 

1 o2 J J 


(IV- 10) 
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(IV-16) 


q = -7 I PjYi 

2 02 j=j J J 


3P: 


with the obvious notation p. = The development of the locally opti- 

mum estimation algorithm will be complete with the evaluation of the 
components In the tg of equation (IV-IO) and the consequenl inver- 

sion to provide 4> . Therefore, proceeding to evaluate the entries in 
(IV-IO), yields 




j=l 


J J 


(IV-l7a) 


■ "T I p^ 

2 aS jii'j ^ 


(IV-l7b) 


<q,q > = <q,q > " I “Pj^ 

12 21 Of^j_|JJ J 


(lV-l7c) 


Note that (lV-l7c) Implies that if p^. has even symmetry and has 

odd symmetry about the mid-point of the summation, then <q q > = 0 ( <y?> 

1 2 J 

has even symmetry if p. does). The conditions required for this to be 
true are: 

1. the antenna selectivity function p(eg) has even symmetry .tbout 
its bores ight, 

2. the derivative of pieg) with respect to the angle off boresight 
has odd symmetry about its boresight, and 

3. +he sampling times are symmetrically distributed about the 
center (boresight) of the stored signal p(0g). 

2 

The only significant problem remaining is the evaluation of <y.> • 

From the defining equation for y,, (IV-14), it is clear that 

J 

<y2> = <[v. 1^1 1 ~\ > - p2 (IV- 1 8) 

J io( ^ J 

where 
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Iflt ) ' ^0 J Ifl' J 2n^ J 


Il( > 


v4 + 


(1V-I^1) 


I ( ) ■ ■ 

Due to the nonlinear nature of - 7 ^, — r , where ( ) •" ( , this inteqr.tl 

has not yet yielded to an exact evaluation. Certainly this integral could 
be evaluated via careful digital computation; however, the value of the 

p4 


integral is a function of the (s igna I -to-noi se) ratio R*: = ^ 

^ J 2o2 2o^ 

= p4 and would, therefore, require on line computation for each sample 
J -2 

value p. and scan estimate R = — - . As this appears to be an unreason- 
ed 25^ 

able computational buroen on the anticipated receiver microprocessor, an.i- 
iytic approximations were investigated. In order to minimise the .jS'.ump I i in-, 
made, power series approximations I 0 Ihe <un< 1 ion wore used IhjI IIk' 

resulting integrals, which were put in the form of the moments of tlio prt li.i- 
bility density function associated with 0 ?(x,y)f *Jre given in terms of 
the confluent hypergeometric function which is expressed either as an 
infinite sum or in terms of exponentials and Ij and Ig. In the first case 
the value of the integral is expressed as an infinite sum of an infinite 
sum, and in the second as an infinite sum of products of exponentials and 
modified Bessel functions. In neither case were closed form solutions 
evident; nor was it clear how reasonable assumptions might simplify the 
expressions. 

Approximations to the values of the integrals can be accomplished 

by examining the nature of the nonlinearity and the parameteri zat ion of 

the Rician density function, i .e. , its functionai character for fiigh and 

low values of r 4, e.g. [7, p. 4|4^ Gaussian for large r 4 and Rayleigh for 
d d 

sma I i Rj (= 0). That this is a reasonable approach is supported by tfio 

fact that the antenna selectivity function Pj is, by design, a highly 

selective function, i.e., its values tend to be relaMvely high in the 

ma inbeam and fall off rapidly to relatively low values in the side lobes, 

e.g.-23db. Thus it seems reasonable to dichotomize the problem into 

these two extremes of Rj and use asymptotic approximations in these two 

cases. 
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Case I: » I 


Bf>caiiso th« final result deuireri is the evaluations of <v^'" • 

J 

given In (IV- 18) which Is the difference in two relatively larq<? numtx'f. 
which may differ only slightly for Rj >> I , care must be taken in iipproxi 
mating the value of the Integral in (IV- 19). For example, assumim 
T,( ) ''''’• 


lo< > 


T “ I for ( ) = 


>> I and assuming the Rician density becomes 


Gaussian with mean P. and variance o* yields <V^> = with <y^> = o^. 

J f J j n j n 

Making the same assumption that = I but retaining the exact form of 

the Rician density yields [7, p. 415] <V?> = + 2o^ with <y2> = 2o^ . a 

J J n 'j n 

variation of 2-to-l from the previous result. The following analysis 
produces an upper bound on the value of the integral which appears to be 
reasonabiy accurate. 

^ function in the 


To avoid squaring the approximation of the 
integrand, use, dropping the subscripts. 


Iq( ) 


[V 


hLl 

Iq( ) 


lo( ) = Cv Ii( ) • 


(IV-20) 


with (IV- 1 9) now becoming 




dV 


(lV-21 ) 


= P ^V I ^ , expectation wrt N 4 (P,o^) [9, p. 42] 


Now introducing the approximation | = l(> — ^) 

Iq^ ' Iq' ' 




N4 


N4 


(IV-22) 


< P(2o2)V^%(2.5),Fi(2.5;2;R2) 
n * * 


(2lt)’ n • 

^ -r- 

4 n 


'^^''^^In(R2/,)[3 ^ + 2R (I + 


vrw 


TqTrW 
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This upper bound Is used In obtaining the upper bound on the ratio 

X = Table (IV-l). This upper bound calculation, in addiii^" 

to several digital computations of the values of the intonral of (IV-h"", 

loads to the choice A - I, I.e., <y^> = for > I . 

— I- j n 

Case II: < I 

.1 

1 1 ( X ) X 

Us I rig the same approach as abfjve, excepf .ipfiroximit I ng "j “ 

for R^ < I , yields 


<v| 44 > 


Io( ) 


20n N4 2q2 


(lV-25) 


p2 - _r2 

S (r-T)2o2e ^ r(3)iFi(3;2:R2) 
20n n 


< 2P2(| + 


Using this upper bound for R^ < I yields 


<y2> a P2(| + r2) 

J J j 


(IV-.’T) 


To use these approximations in (IV- 1 7) assume that there are — samples 
taken In each of the to-and fro-scans (K/scan) and that Rj < I for I < j < f 
and > I for t+\ < j < K. Therefore 


j=i J J J j=ii J " 


(TV-?^J) 


I 


K 


j=| ■' J I I 


(IV-2bb) 




and for Rj << I for I < j < ^ 


T *2 


<qi> = 2 R 2 V 

maXj^^l J 


(IV-?bc) 
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Tablo (IV-I) 




4 

9 

16 

25 

% 

49 

64 

81 

100 


UPPCR BOUND for <y?>/oj^^ for R-“* ;; I 

I .2/411 
1.44897 
1.47941 
I .48853 
1.49248 
1.49477 
I . 49609 
I .49695 
1.49753 
1.49792 
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There Is some computational evidence fhat <0|> - 60 R^ ‘»n.i !uf ffu>r 

” max 

investigation of the validity of this simple approx i m.i1 ii*n will he e ' we.*. 
Conflnulng with the evaluation of the elements of ♦ leads to 


<qi> • -V t I ♦ R^) + y p^o2} 
On j = | J J J j=h| J " 

* -T ^ 5^ + I Pj} 

On j=| J J J j=Z+l 


(IV- 20a) 


(IV-20P) 


(IV-26C) 


(IV-2bd) 


and for the present study p. and p. will he assumed to have fhe requirr>c1 

J J 

symmetry to have <qiq 2 > = 0. None of the forms of tjquafions (IV-25) or 
(IV-26) present significant computational dlfflculfy and the sensiMvity 
of system performance to the accuracy of approximating those terms will 
be Investigated via simulations. Returning to the estimate ^ of (IV-l) 



(IV-?/a) 


(IV-27h) 


and the Individual error estimates are 


*8 = <qf> ' <,'}> • ^ itfr- '’j 


(IV-PH.J) 
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Il( ) 


e = 9 Q? ” o • -^9 y p . {V , ■ ' : "• P , } 

a <q^> <q2> ^,2 ^A^Pj j Iq( ) 


(IV-28b) 


where the indicates parameter evaluated at the estimate from the 

previous scan, e.g. , P. = oP. = apC0,(T.) - 0(t, )H. The above error 

J j Aj k 

may be used to update the tracking receiver estimates of the angular 
position of the A/C, 0, and the peak amplitude a of the direct path signal. 

Note that if the sums over samples for which Pj < I, i.e., I < j < t, 

are neglected In computing <y^>, the resulting equations for e and e are 

j 0 o 


K 

I Pj{V 

j = ! ^ 


IiLl 

j Iq( ) 




(IV-?9a) 


K T / ^ / K 

' = I 6.(V. P.) / I 

“ j=| J J Io< > J ' J=£+| J 


(lV-29b) 


where the denominator sums are nearly cons*'ant, I.e., t constant, over a 
fairly wide range of and thus may be precomputed if on-line computa- 

tional power is limited. 

Preliminary results of simulations are presented in Figure IV- I . It 
should be noted that, at this point, no attempt has been made to assess 
the accuracy of the approximations in <qf> and <Q 2 > or fhe effects of these 
inaccuracies on the performance of the estimation algorithms (sensitivity 
analysis). Again it should be pointed out that the error statistics (one 
o values) given here are for the raw error which contains none of the 
smoothing anticipated when a tracking algorithm is in use e.g. those 
suggested in Section II of this report. 


As was pointed out in Sections II and III, if e = 0 the residual 
mean square error associated with the estimate e is <J>o, i.e., 

0 


<ee^> 


e=0 



(IV- 30) 


<q2> 
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(RM5 0 ERROR IN DEGREES) 



I 


Using the simple approximation following 


and, therefore, <e^> 
(lV-25) yields 



<e^> = 


60R2 

max 


or 


(lV-3la) 


-oh 0.13 
<e^> = 


max 


(IV-5lb) 


for theoretical root-mean-square values of 0.013 and 0.051 for of 100 

max 

(20db) and 6.3 (8db) respectively. These numbers are very close to those 
obtained in simulation (0.012 and 0.047, see Figure IV-I) and tend to 
show agreement between these theoretical developments and computer simu- 
lations. However, as the simulation results appear to produce smal ler 
root-mean-square errors than the theory would indicate possible, there 
appears to be some need for further refinements of this analysis. Addi- 
tionally, the simulations produce a significant bias in e which Is, 

n 

of course, counter to the theoretical development of an unbiased estimator, 
i.e., <e > = 0. There is some evidence that this bias is due to a very 

O 

remote characteristic of the specific computer program being used to 
generate the random noise in the simulation. This problem is now under 
study and alternate noise generating programs are being considered. 

Comparison of l-F & Envelope Detectors 

From Section III Equation (III- 14), the mean square error of the i-f 
detector is given by 


<®^|F 


2No 


e=0 


7' p'i:e,(T)]dT 


(IV- 32) 


which can be closely approximated by 
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<©2 > — 
IF 


2N, 


a^AT I 

j=l J 


2NqB 

«2 


K 

J=l J 


A 9 r\ 

= I P]) 

Op j=i J 


-1 K 

= (2r 2 y p2) 

max 


-1 


(IV-35b) 


where At = Is the reciprocal of the i-f bandwidth and 2NoB " i-otn- 

u n 

paring this with the results for the envelope processor, equation (IV-2bc), 


^ A 

<e^> 

0 E 



= (2R2 

max 


K 




(IV-34) 


and the ratio for the two cases becomes 


<;2> 


0 IF 


/ <«0^E = 


K 

I 

j=£+i 




K 

^Pj - 
j = l ^ 


(IV-35) 


For large signal-to-noise ratios, i.e., 


K 

I 


Ip^ 

j = l 


j=£+i 


mean-square performance of the two algorithms become nearly equal, 
r 2 = 100 (20db) the ratio is 0.94. Also the numerator of (IV-35) 


and the 

e.g., for 
is a 


max 

lower bound on the actual value, and yet it can be shown by actual computa- 
tions that the value of the .ratio of mean-square errors is bounded above 

by unity for R^ sufficiently large, e.g., R?_ > 6.3 (8db). As the 

max max 

performance of the two processors is very nearly the same for reasonably 
high signal-to-noise ratios, all future efforts will be concentrated on 
the envelope processor which is much simpler to implement. One plausible 
explanation for this perhaps surprising result is that there is no informa- 
tion concerning the A/C angle in the phase of the i-f signal as modeled. 

A brief elementary study showed that, in the far field of the antenna, this 
signal model was appropriate. 
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Plans for Future Work 


Work will continue in many of the areas discussoil in this •;c'c*i.''' 

with respect to the envelope detector processor, e.q. ().ir.ime1 r i* s»'ii- i- 

tivity study via simulations, further analytic ref i nornent of the houniJs 

on the formulation of <y^>, etc. New work will include analytic and 

J 

simulation effort in the following areas: 

1. Multipath suppression including use of both random process .nul 
unknown deterministic signal models of the multipath, 

2 . Revise signal model and resulting processor derivation to include 
logarithmic l-f as well as alternate types of detectors, e.q., 
linear/logarithmic squared amplitude envelope, 

3. Investigate the feasibility of making the envelope processor 
adaptive to the antenna selectivity function; also look for 
robust designs which might be insensitive to this possibly vari- 
able signal feature without significant loss of error performance, 
e.g., the square gate receiver used in our simulation. 

Most of this work is scheduled for completion prior to the start 
of programming the microprocessor, i .e. , June I976. All of fhis porNon 
of the project should be completed by mid-July I97fi. 
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V. CURRENT INVESTIGATIONS, SYSTEM CONS I DERAT IONS AND RECOMMENDAT I 


All effort currently is related to optimal onv('lope process i no . 
Theoretical investigations have been discussed and arc‘ summarized lx ‘lew; 

1. Multipath-adaptive processing of the received signal envelop<!. 

2. Comparative evaluation of finte-word length digital processors 
designed to accept linear envelope, log envelope and squared amplifudt' 
envelope signals. 

3. Comparative evaluation of the layered tracking algorithm dis- 
cussed with structures of the ful ly-recursi ve design. 

4. Tentative consideration of algorithmic requirements for a 
beamw i dth-adapt i ve feature. 

Simulation work will follow the theoretical studies listed; in addition 
the receiver evaluation filters will be added to the simulation to l<jcili- 
tate evaluation and comparison of results with those of other studitjs. 

System design of the prototype signal processor to bo I I ighl-teslod 
has begun. Microcomputing equipment has been ordered, though previously 
unanticipated delivery delays has forced some change in the project 
schedule; a revised schedule is given in this report. A longer period 
of more extensive pre-flight exercising of the system than originally 
scheduled is desired, and it appears this may not be possible in the 
current funding period. This is discussed further below. 

In general the software in the microcomputer must accomplish foui 
distinct tasks: 

1. Input conditioning and storage 

2. Algorithmic calculations 

3. Output updating and posting 

4. Allocation of the machine resources to the above 3 tasks (i.v.. 
the executive program). 

The input conditioning and storage functions will be served by a diro< t- 
niemory-access (DMA) controller, which will autonomously saripl.' Hir iriuul 
analog (envelope) signal over the proper intervals ( dotermi ned tiy Itn- 
prediction 6(k|k-l) during the TO and ERO scans, perform ana I oq- to-d i g i l.i I 
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conversion of these samples and then sfore the data in prescribed block-, 
of memorv (RAM) in thn machine. An interrupt from thi> 0M/\ cent r.’* I lot 
followifiq storaqe ol data tor the T(' and FRO scan*. euMi will eiiahli' 
processing of the now data (when otherwise timely). 

The algorithmic calculations use only memory tor source and sync 
ot data and hence will run in the "background" under executive control 
(of initiation). 

the output estimcjte is obtained by "rrvjrking" 1h(> ext rapo la t i n«i 

estimate at the occurrence ot the Barker correl.Oion peak (t pulse) 

o 

that indicates the arrival ot the next scan. A default trigger will 
provide this function (with suitable advisement) in the event of loss 
of sync signal. The angular coordinate estimate is then passed both to a 
serial output port tor external distribution and to a parallel latfhod 
output port tor local digital display. 

The executive program takes all external timing and channel identi- 
fication information and performs the scheduling function in the miichino, 
inciuding tor example; 

1. "Mark" the extrapolating estimate and terminate extrapo I at i or. 
when the new sync pulse arrives (it sync signal is lost, mark-by-defau I i 
the extrapolating estimate but do not terminate extrapolation) . Outpuf 
the "marked" value of the angular coordinate estimate. 

2. Adapt the algorithm to the elevation or azimuth channel, as 
appropriate to the upcoming scan; 

3. Initiate processing ot TO scan data after storage is complet<'. 

4. Initiate processing of FRO scan data and subsequent estimafr- 
extrapol ation after storage ot FRO scan data is complete. 

5. Resolve ambiguities in timing and syncing associated with turn-r.n, 
signal loss and recovery, etc. 

6. "Idle" when all scheduled calculations are complete. 

if the input data rate (multichannel) is too high, the executive prooram 
should also selectively edit the input data stream in an acceptable manner 
to not exceed the processor thru-put rate. 


31 



Above, the advisability of a somewhat fuller lab testing and I'le- 
Hight exercising of the prototy,jo was mentioned (in .omii'cMv'M ,\ I ■ 
computer delivery del.iys and associated schc>dule imp.ictst. Ilii ;’i i i. d 
njason for this is ttuit 1tie prototype system, <is it is evolving, will 
involve many different functions, only a few of wtiirh (I tie a Igoi i | luni. 

' a I cu I a t ions ) will if have been feasible to tost in ttie large FORTRAN 
simulation. It would seem prudent, if possible, fo Inst a I I functions 
prior to interfacing the prototype with the Phase? Ill Receiver. Relevant 
als(? is the available electrical and functional output of the Ph.ise III 
Receiver, its full definition, and the potential necessity to install 
line drivers in the Phase III Receiver enclosure to send signals by 
coaxial cable to the prototype system. Our recommendations for future 
work are essentially that provisions be made for those tasks, particularly 
the system test. A real-time simulation of the Phase 111 Receiver can 
bn nconomically doveloped in the mi vti •ctirii|Mj I lu d' -v 1 1 n t sy in Iti'' 

l.ib and used to exercise the complofod profolvpn sysfnm in tofo tlirougli 
it - designed hardware interface. This would h(?ln to insure a smooth 
iti+ngratlon with thc' Phase III Receiver .ind to ciuar.in lt?e the expected 
performance at flight test time. A proposal for project continuaticr 
along these lines will be submitted. 
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